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Abstract 

We investigate multichromophoric energy transfer allowing for bath-induced fluctuations at 
different sites to be correlated. As a prototype system we consider a light-harvesting antenna 
surrounding a reaction center. We show that the interplay between quantum coherence and 
correlated fluctuations can generate a room temperature transfer process featuring a marked 
dependence on the degree of symmetry and derealization of the initial exciton state. Our work 
illustrates how these quantum features could support fine tuning of energy transfer efficiencies 
in closely-packed natural and artificial light-harvesting complexes. 
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Remarkable experimental advances in multichromophoric assemblies have raised the fascinat- 
ing possibility that quantum coherent dynamics may be relevant in photo synthetic energy trans- 
fer, even at room temperature.^^ Key to the survival of quantum coherence in this temperature 
regime seems to be the emergence of correlated energetic fluctuations between different chro- 
mophores. 6 This could be due, for instance, to closely spaced pigments sharing the same envi- 
ronmental modes. 2 3 There is thus a growing need for theoretical descriptions of energy transfer 
beyond the common assumption of independent environments for each pigment.™ Indeed, a num- 
ber of studies have addressed this issueJ^HHl though the subtle interplay between coherence and 
correlated fluctuations in influencing transfer efficiencies remains largely unexplored. 

The prospect of environment-protected, room-temperature quantum coherence in photosyn- 
thetic systems has revived the longstanding question on the functional role such quantum phe- 
nomena may play. In particular, it is natural to ask whether the interplay between coherence and 
correlated fluctuations provides the system with an efficiency control mechanism that it could not 
otherwise achieve.^ In this work, we address this question with the aid of a model of multichro- 
mophoric energy transfer under the influence of spatially correlated energetic fluctuations. We 
show that a fundamental signature of the interplay between coherence and correlations is a quan- 
tum yield with a conspicuous dependence on the phase information embedded in the initial exciton 
state. In particular, the energy transfer efficiency can distinguish the degree of symmetry or asym- 
metry of the initial exciton state, as well as its derealization length. Such a dependence produces 
an ordering of the excitonic eigenstates from higher to lower efficiency as a function of increasing 
energy, a feature which does not exist under local environments. Hence, coherence and correla- 
tions can indeed equip large, densely packed photo synthetic systems with a unique mechanism to 
control their transfer efficiency. 

As a prototype system we consider a light-harvesting 1 (LH1) antenna connected to a reaction 
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Figure 1: Schematic of the LH1-RC core of purple bacteria based on Ref.^ Arrows indicate 
the transition dipole moments for each bacteriochorophyll (Bchl), as given in Ref.^ The special 
acceptor pair in the center (red dots) defines a set of preferred axes of symmetry for the electronic 
couplings between the center and the LH1 Bchls.^ 



center (RC) as in purple bacteria™ (see|[figure][l][]l[). We focus on this unit as it exhibits key 



features, such as symmetry in its geometry and interactions, representative of LH units present in a 
large variety of natural species 21 and artificial systems;^ such symmetries are expected to play an 
important role in the energy transfer process. Evidence of coherence in the LH1-RC core is sparse, 
though its excitonic states can be delocalized over a large number of pigments at low tempera- 
tures. 23 Moreover, it has been shown that correlated energy fluctuations among several Bchls in the 
LH1 ring must be included to accurately reproduce spectroscopic properties, 24 pointing again to 
the importance of environmental correlations in these densely-packed pigment-protein aggregates. 
The Hamiltonian describing the electronic excited states of the M strongly interacting pigments of 
the LH1-RC core is given by H s = £m=i (E m + 8E m )a+a m + Zn> m Vmn(o+a~ + a+a m ), with 
h = 1 . The operator <7+ creates a localized excitation at site m, with an energy separated into an en- 
semble average E m , and a deviation 8E m that accounts for diagonal static disorder (see supporting 
information). On-site energies E m and electronic couplings V mn have been derived from RefsM^-^ 
and they reproduce the observed exciton spectra.^ In contrast to the common assumption of lo- 
cal, independent environmental modes in contact with each pigment, we consider here a single 
thermal bath described by H B = £ k (O^byb^, where (Z?k) are the creation (annihilation) opera- 
tors of bosonic modes of frequency 0)^. Each pigment is assumed to have a position dependent 



coupling to the bath of the form (g m )k = |(gm)kk' k ' r " ! , with r m the position of site m with respect 
to each phonon mode. 26 This position dependence will give rise to correlated energy fluctuations. 
For simplicity, we consider isotropic and identical coupling strengths for all sites, |(g m )k| = |gk|> 
with associated single-site spectral density J? {(Q) = Ek |gk| 2 <5(co — (Ok), and reorganization energy 
Er = JqcIco J?{(q)/(Q, quantifying the form and strength of the system-bath coupling, respectively. 
Hence, the system-bath Hamiltonian reads H$b = Lm=l °m°m ®B m (r m ), with bath operators writ- 
ten in the interaction picture as B m (r m ,t) = £ k |g k | ^ k e'^ t+k ' r '^ + b it e^ i ^ 01kt+k ' r ' n ^ . 

Here we consider the limit of weak system-bath interaction for which the energy transfer dy- 
namics can be described with a master equation derived within the Born-Markov and rotating- 
wave approximations:^^ p = — i[H$ + His,p] + 2$(p). Here, p is the reduced system den- 
sity operator while His is me Lamb-shift, which will be neglected. The dissipator, 3>{p) = 
HcoHm,nYmn(o) : d m n)^mn(<X),p), accounts for the effects of both local and correlated fluctuations 
on the system. The operator ^ mn ((0,p) = A n ((o)p(t)Aj n ((o) — \{A^ m ((Q)A n ((Q),p}, is written in 
terms of Lindblad operators A m (co), with {00} the frequency spectrum given by the energy differ- 
ences between single-excitation eigenstates of Hs- The rates y m n{o),d mn ) = f^ oa dte m ' t C mn {t ,d mn ) , 
where d mn = \r m — r M | is the separation of sites m and n, are defined in terms of bath correla- 
tion functions C mn {t,d mn ) = (B m (r mi t)B n (r nj 0)} B , with the average taken over the thermal equi- 
librium state of the bath. Evaluating the correlation functions, we obtain a formula relating 
the spatially-correlated dephasing rates to the on-site rates, Ymn(co,d mn ) = ^(k(a))d mn )Ym?n(o)), 
where the single-site rate is y mm (co) = y((o) = 2n ( /{\(o\)\N{— m)\, with N(co) being the ther- 
mal occupation number. The function 3f (z) depends on the environment dimension, exhibiting 
similar decaying behaviors in two and three dimensions. 26 Here, we assume a two-dimensional 
bath such that (z) = Jq(z), with Jq(z) being a Bessel function of the first kind. The func- 
tion k(co) defines an effective bath correlation length, which as a first approximation we assume 
to have a frequency-independent form k(co) = 1/Rb, accounting for a phenomenological cor- 
relation length Rg. 9 We can now simplify the dissipator as 3>{p) = L ffl y(fi)) [Lm^mm(c0,p) + 
Y. n ^mJo(dmn/RB)^mn((0 1 p)], where the first term accounts for independent, on-site fluctuations, 
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and the second for correlated processes. Using the property Y^mn^mniS^iP) = 0, we can demon- 
strate that for Rb — > °°, we have 2$(p) — > 0. In this limit all sites are fully correlated, and the 
single-excitation subspace of H$ is effectively decoupled from the bath. In the opposite limit 
of vanishing Rb, we recover the case that each site experiences only local fluctuations, @(p) = 

r(co)^ mm (co,p). 

Over the course of the system's evolution any initial excitation will eventually be either dissi- 
pated or trapped, with respective rates Y m and K m for each site. These two incoherent processes 
can be accounted for by the inclusion of a non-hermitian operator, i/ pro j = — *X m (r m + K m )o^o~, 
into the system Hamiltonian.^^SI efficiency of transfer (quantum yield) is then defined as the 
total probability that the excitation is trapped, r\ = Y.mJo K m (m\p(t)\m)dt, where \m) = c+|0). 
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Figure 2: Transfer efficiency at room temperature (293 K): (a) versus on-site reorganization energy 
En in the absence of correlated dephasing. Er increases from right to left along the horizontal axis 
to allow easier comparison with below; (b) versus Rb for Er = 100 cm -1 . 

Results - In the following, we assume that all sites dissipate at a rate T = 1 ns -1 ,^ with the 
exception of the special pair at the RC, which traps excitation at a rate K = 4 ps~ 1 . 2 ^ For LH com- 
plexes the system-bath interaction is usually characterized by an Ohmic spectral density with either 
Drude^or exponential cutoff;^ we assume a Drude form, ^ (fi)) = 2Er(d c (q/ n((0 2 + CO 2 ). Typical 
Er values are estimated to be between 100 cm - and 200 cm 1 31 for LH1, while the bath cutoff 
frequencies (O c used to reproduce experimental data range from 50 cm^^to 1000 cm -1 . 24 Here, 
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to be consistent with our weak system-bath coupling model and to unambiguously illustrate how 
initial state distinguishability can arise through spatial correlations, we consider a fast-modulated 
bath with (O c = 300 cm -1 , and a maximum Er = 100 cm -1 . The ring sites are numbered from 
1 to 32 (see [figure][l][]l[ ) and we assume pure initial states, having no overlap with the RC, of 
the form = (Vv^)£yLi(±l)' / !./'}> where 1 < m < 32. The initial state is thus symmetri- 
cally or asymmetrically delocalized over the ring. In particular, we focus on four representative 
states l^g), |*Pg ), 1^32), and ^32), which in the case of pure coherent evolution (Rr — > 00 or 
Er = 0) have, in general, different transfer efficiencies. 18 Further, | v J / 32 ) (1^32)) has a large degree 
of overlap with the lowest (highest) eigenstate of Hg. Since the dynamics of p varies with respect 
to pigment-exchange, the results presented correspond to an average over all initial states with 
derealization over m consecutive sites. 

In [figure] [2] []2[ a), we plot the transfer efficiency, rj, for each initial state as a function of the 
system-bath coupling quantified by Er, and accounting only for local dephasing, i.e. y mn (o)) = 
for m 7^ n. For small Er, r\ varies significantly between the different initial states, consistent with 
the coherent limit.^- 8 -' As Er increases, exciton relaxation and dephasing take place and above Er > 
10 cm -1 (dashed line), the efficiency becomes practically identical for all initial states considered, 
indicating that energy transfer is dominated by incoherent processes. Notice that Er ~ 10 cm -1 
is comparable to the electronic interaction between a Bchl in the ring and the special pair in the 



center. 18 The situation differs for correlated dephasing, as shown in [figure] [2] [] 21)). Here we fix 
Er = 100 cm -1 , which falls in the ranges of reported values for an LH 1,^3 and plot rj as a function 
of the bath correlation length Rr. For large enough values of Rr the efficiency distinguishes again 
the initial state meaning that coherent dynamics is at play. Two key features can be drawn out from 



the comparison of [figure] [2] [] 2 (a) and (b). First, for moderate values of Rr the transfer efficiency 
clearly separates symmetric from asymmetric states - we shall shortly discuss the role of coherence 
and correlations in such distinguishability. Second, for an extended system, increasing Rr does not 
generally correspond simply to a renormalized coupling to an effective local environment at each 
site. This contrasts with the intuition to be gained by considering the two-site problem (M = 2). 
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In this case, we again use Y*mn^mn{P-iP) = to rewrite the dissipator as S>(p) = L£o7( ft) )[l — 
Jo{d / RB)\Y,m=\ 2^mm{®iP)i where d is the inter-site separation. Here, a finite Rg is equivalent 
to considering each site to be independently coupled to its own local bath, but with an effective 
reorganization energy E R = E R [l — J^{d / Rg)]. 




Figure 3: Main: (Left axis) Room temperature transfer efficiency for |*Fg ) against Rg, with Er — 
100 cm -1 : (•) assuming local dephasing with an effective Er = Er[1 — (J(d mn /RB))}', (o) including 
only correlations satisfying J(x) > 0.7; (line) including all correlations. (Right axis) (+) Effective 
number of correlated neighbors to each site. Inset: Tj versus Rg for the state l^i). 
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Figure 4: (a) Efficiency for asymmetric and symmetric states as a function of derealization length 
m for Rg = 40 A. (b-c) show r\ in the absence of correlations (Rg = 0) for asymmetric (b) and 
symmetric (c) states, (d) Efficiency for eigenstates of H$ having vanishing overlap with the RC for 
Rg = AQA (o), andi? s = (o). In all cases E R = 100 cm" 1 . 

Interpretation of the effects of correlated dephasing therefore requires us to go beyond the two- 
site model. We analyze three different regimes: (i) Rg < d m i n , with d m i„ ~ 10 A the minimum pair 
distance in the system; (ii) Rg 3> d max = 2R, where R ~ 50 A is the radius of the ring; and (iii) the 
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most interesting, intermediate regime d m i n < Rg < d max . In the first two cases the two-site intuition 
can be applied as fluctuations are only very weakly correlated in case (i), while they are almost 
completely correlated at all sites in case (ii). In neither case is there significant variation from the 
average correlation over all m 7^ n pairs, (Jo(d mn /RB))- Hence, the dissipator has approximately 
the same form as the two-site case with an effective reorganization energy Er[\ — {JQ{d mn / Rb))\. 
This simplified dissipator does not reproduce the transfer efficiency in the intermediate regime for 
delocalized states, as can be seen in 



[figure] [3] [] 3 where we plot r\ for the initial state |*P 8 ) as 



an example. Interestingly, for the state l^i), which has no quantum superpositions, the effective 



local-bath approximation captures the efficiency behavior across the whole range of Rg (see [fig- 



ure] [3] []3 inset). The key feature, then, is that there is a crucial interplay between coherence, bath 



correlations, and initial state properties in defining the transfer efficiency. To understand this in- 
terplay we can follow an analysis within the quantum jump picture. The system dynamics can be 
seen as a series of coherent evolutions interrupted by incoherent jumps between single-excitation 
eigenstates, described by the action of the operator Y mn {(o)A m {(o)p{t)AJ n {(o). In the presence of 
bath correlations, the overlap with the trapping sites due to such incoherent jumps carries informa- 
tion on the quantum phase differences between sites m and n, through a factor (±l) m+ '\ while in 
the absence of correlations it does not. Hence, after a jump due to a correlated bath, the subsequent 
coherent evolution is weighted by a phase factor that is always 1 for symmetric states, but can be 



1 or —1 for asymmetric states. This gives rise to the symmetry distinguishability observed in [fig- 



ure] [2][]2"b). In short, correlated fluctuations promote an initial state symmetry dependent energy 



transfer towards the RC. 

Further insight into the behavior of 77 in the intermediate regime (iii) can be gained by asso- 
ciating Rb to an effective number of neighbors correlated to each site. We estimate this number 
by neglecting all correlations below a cutoff y (Jo(d a /RB) < y), chosen to be the largest value that 



still gives a good agreement with the exact results. For the relevant points in [figure] [3] [] 3 we take 
y = 0.7, which corresponds to neglecting all correlations between pairs whose distances satisfy 
d m n > 1 In Ref. 24 correlations among 5 consecutive sites were included to fit measurements 
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of optical properties of LH1 samples. From the above analysis, we can associate this number to a 



correlation length of Rb ^ 20 A (see [figure] [3] []3 ). Notice that at this moderate correlation length, 



r\ already starts distinguishing states according to their symmetry (inset in [figure] [2] [] 2 b)). 

With Rb now associated to an effective number of correlated neighbors we are led to conjecture 
that, for a given Rb, the interplay with coherence will be such that r\ should increase as a function 
of the derealization length (m) for asymmetric states, while it should decrease for symmetric 
states. Further, such behavior should happen only up to an optimal delocalization length m c , 
beyond which the efficiency should not considerably change. This is indeed the behavior seen in 
[figure] [4][]4[ a), where we plot r\ for asymmetric and symmetric initial states as a function of m, 



for an intermediate value Rb = 40 A with associated m c ~ 12 (dashed line). These results confirm 
that the interplay of coherence and correlations allows the system to exploit not just asymmetry, 
but also exciton delocalization. We find, however, that the optimal m c for which the efficiency 
plateaus increases only slightly with larger Rb, and is never greater than m c ~ 16. Hence, m c 
is not simply dictated by the length of correlations, but is mainly defined by the symmetry of 
the electronic interactions in the LH1-RC core. Further, the marked splitting between the two 
sets of states around the efficiency associated to m = 1 clearly illustrates the symmetry-dependent 
energy transfer discussed above. To highlight the difference to the uncorrelated environments case, 
[figure] [4][]4[ b-c) show rj as a function of m when Rb = 0. Importantly, although the variation in 



r] is negligible, its functional form for each set of states is the same as that observed under purely 
coherent evolution.^ 

In the actual light-harvesting process the implications of the above discussed phenomena may 
be rather subtle. To illustrate this, in [figure] [4][]4[ d) we plot rj against energy for all eigenstates 
having vanishing overlap with the RC. The presence of correlations is manifested in an unexpected 
ordering of the eigenstates whereby r\ decreases monotonically with increasing energy. Such order 
is not observed when correlations vanish. In a functioning photo synthetic system the LH1 donor 
state will be a thermal distribution about the lowest energy levels.^ Our results indicate that, 
without correlations, such a distribution will be as efficient as one about higher states. In contrast, 
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if correlations are present, a thermal distribution about low-lying states would be the best option. 
This suggests that the interplay of coherence and correlations could indeed help the natural system 
to differentiate the best energy transfer pathway. 

In summary, we have shown how coherence and bath correlations could equip light-harvesting 
systems with a mechanism to exploit initial state phase information at room temperature. Although 
our conclusions are based on a simple model, we conjecture that initial state distinguishability in 
energy transfer may be a general feature when spatial and/or temporal environmental correlations 
are present. Hence, experimental explorations of such phenomena could help in elucidating both 
microscopic details of the environmental correlations as well as the advantages of coherent dynam- 
ics in natural and artificial light-harvesting systemsJ^EH 
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